Probing the tidal disruption flares of massive black holes with high-energy neutrinos 
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The recently discovered high-energy transient Swift J164449. 3+573451 (Sw J1644+57) is thought 
to arise from the tidal disruption of a passing star by a dormant massive black hole. Modeling 
of the broadband emission suggests the presence of a powerful relativistic jet, which contributes 
dominantly to the observed X-ray emission. Here we suggest that protons can be accelerated to 
ultra-high energies by internal shocks occurring in the jets, but their flux is insufficient to account 
for the observed flux of ultra-high energy cosmic rays. High energy protons can produce ~ 0.1- 
10 PeV neutrinos through photomeson interactions with X-ray photons. The large X-ray fluence 
(7x 10 _4 ergcm -2 ) and high photopion efficiency, together with the insignificant cooling of secondary 
mesons, result in bright neutrino emission expected from Sw J1644+57 if the jet composition is 
matter-dominated. One to several neutrinos may be detected by a Km 3 -scale detector from one tidal 
disruption event similar to Sw J1644+57, thereby providing a powerful probe of the composition of 
the jets. 

PACS numbers: 95.85Ry, 98.70Qy, 98.70Sa 
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Massive black holes are believed to reside at the cen- 
ters of most galaxies and the vast majority of them are 
considered to be dormant. It was long predicted that, 
if a star of mass M* and radius i?* passes occasionally 
within the disruption radius Tt — ^(-Mbh/M*) 1 ' 3 of 
the dormant black hole (where Mbh is the black hole 
mass), the star will be torn apart by gravitational tidal 
forces, leading to a transient accretion disk and a bright 
panchromatic flare [lj. Candidates of such tidal disrup- 
tion flares (TDFs) have been suggested Relativistic 
jets are expected to form in such accretion disk system 
and may produce observational phenomena as well[3j. 

The high-energy transient Swift J164449.3+573451 
(hereafter Sw J1644+57) was discovered by the Swift 
satellite at 2011 March 28 12:57:45 UT0|. The facts that 
it occurred near the center of a galaxy, no archival X-ray 
emission before detection and its analogy with blazars 
in the X-ray emission suggest that it is mostly likely to 
be a TDF event 0-0). At redshift z = 0.3534 (corre- 
sponding to a luminosity distance of d^ — 1.8Gpc)@, 
the early X-ray flare emission reached a peak luminos- 
ity of Lx — 4.3 x 10 48 ergs _1 (isotropic equivalent) and 
then transits to a low state with a median luminos- 



tic jet of a bulk Lorentz factor of T ~ 10 pointed towards 
us 0, [E| . Modeling of the broadband spectral energy dis- 
tribution (SED) also requires a powerful relativistic jet 
with T ~ 10, which produce dominantly the observed X- 
ray emission^, @. A relativistic jet is also required to 
explain the radio transient [![. However, how the jet is 
launched is not well-understood. It is believed that the 
composition of the jet, whether it is matter-dominated or 
magnetic field-dominated, is crucial to unveiling the for- 
mation mechanism of the jets. In this paper, we suggest 
that jets in TDFs can produce bright emission in high- 
energy neutrinos if the jet is matter-dominated, and thus 
neutrino observations provide an important tool to diag- 
nose the jet composition. 

Internal shocks and particle acceleration. — We con- 
sider that a TDF event produces a relativistic matter- 
dominated wind of luminosity L w ~ 10 49 ergs _1 , moving 
with a bulk Lorentz factor T ~ 10. Variability of the 
source on time scale t v , resulting in fluctuations in the 
wind bulk Lorentz factor F on similar time scale, would 
lead to semi-relativistic internal shocks[9[ in the ejecta at 
a radius 



2.96 x lO^ergs" 1 in 0.4-13.5 keV over a 



10 6 s Correcting for the live-time frac- which is wel1 above the photosphere radius at R ph 



ity of Lx 
time AT < 

tion of the observation, the total unabsorbed fluence is 
S x = 7.1 x 10- 4 ergcm" 2 in the observed 0.3-10.0 keV 
band[J]. This gives an estimate of the total isotropic 
equivalent energy of Ex — 3 x 10 53 erg in the 0.4-13.5 
keV rest frame energy band. The observed minimum 
X-ray variability time of Sw J1644+57, t v ~ 100 s0, 
constrains the size of the black hole under the assump- 
tion that the size of the central engine determines the 
shortest variability and suggest an upper limit on the 
massive black hole mass Mbh < 8 x 1O 6 M . The ob- 
served peak luminosity is super-Eddington and requires 
a strongly anisotropic radiation pattern with a relativis- 



R ~ 2T 2 ct v = 6 x 10 14 r 2 ^ !2 cm, 



(1) 



aTL w /(4TrT 3 m p c 3 ) = 1.2 x 10 13 L Wi 4qT^ cm, where ctt is 
the Thomson cross section. We use c.g.s. units and the 
denotation Q = 10 X Q X throughout the paper. 

Denoting es as the fraction of the wind kinetic 
energy converted into magnetic fields, we have 
a magnetic field B' = {%TTe B L w /Ai:R 2 T 2 c) 1 / 2 = 

1.3 x 10 3 e^ /2 _ 1 i^ // 4 9 i?r 4 1 8 r^ 1 G, where the prime symbol 
represents quantities measured in the comoving frame of 
the shock. It is assumed that internal shocks accelerate 
protons with a spectrum dn/de p ~ e~ 2 , where e p is the 
proton energy in the observer frame. The maximum 
proton energy is set by comparing the acceleration time 
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4;c = as p /(eTB'c) = 860a( T ^ v )e B ]L 2 1 L w 1 {g R 14 . 8 s 
with the shock dynamic time = R/Tc = lO^ii^s, 
where a ~ a few, describing the ratio between the acceler- 
ation time and Larmor time. This gives a maximum pro- 
ton energy of e ma x,dyn = 2.4 x 10 20 a- 1 e 1 ^ 2 _ x L 4 J 2 4g T ^ 1 eV . 
The maximum energy is also limited by the 
cooling time of protons. The synchrotron cool- 
ing time is t' syn = 6irmpC 3 T / ' {a T m 2 e e p B 2 ) = 
240e s ;_ 1 L- 1 49 i? 14 . 8 r3( T ^)- 1 s, which gives 

a maximum proton energy e m ax,syn = 0.5 x 

lO^a^^e^L^lgtl^Tl^eV. Another process 
that may prohibit the acceleration of protons to ultra- 
high (UHE) energies is the photopion cooling loss. 
UHE protons of energy e p interact with soft photons 
with energy e 7 = 0.15GeV 2 r 2 /e p = 0.l5T 2 e~ \ eV, 
which locate at near infrared (NIR) band for typical 
jet parameters. The number density of NIR photons 
in the comoving frame is n^ IR = Lnir/(47t i? 2 Tce 7 ) = 
2.5 x 10 14 L N i R:44 rr 5 ^2(e 7 /0.15eV)- 1 cm- 3 , where 
Lnir — 10 44 ergs~ 4 is the NIR luminosity at times 
a few days after the initial trigger 4]. The photo- 
pion cooling time in the comoving frame is t' vl = 



energy injection rate in protons is 



r?e 2 (e 7 /0.15eV)s, 



1/ \a PJ n' mn cKp^) = 1200Z/ NIR44 
where K pl is the inelasticity and <j pl = 5 x 10 _28 cm -2 
is the peak cross section at the A resonance. By 
equating t' acc with t' , we get the maximum proton 
energy limited by the photopion cooling process, 
f - 1 S x ^() 20 ^v- 1 f 1/4 T 1/4 T~ 1/2 r 5/ V /2 eV 

Thus, internal shocks in TDFs can accelerate protons 
to energies above 10 19 eV, in support of the earlier 
suggestion that TDFs can produce ultra-high energy 
cosmic rays (UHECRs) [10| . However, the flux of such 
UHE protons contributed by TDFs in the universe is 
insufficient to explain the observed flux of UHECRs, as 
we show below. 

UHECR flux— The Swift BAT, with a field of view of 
Sfov — 47r/7 sr, has detected one such event in a time 
T ~ 7years at a flux that would have been detectable up 
to a luminosity distance of d max — 5Gpc [4[. Therefore 
we will assume that the rate of Sw J1644+57-like event 
is 



R 



4tt 1 



1 



S F ov T (4/3)^ 



2 x 10 _12 Mpc~V _1 (2) 



and the energy injection rate in X-rays is ix = REx = 
6 x 10 41 ergMpc _3 yr _1 . Following ref.[ll|, the total en- 
ergy in accelerated protons E p can be parameterized by 
E p = £ p Ex, where the non-thermal baryon loading fac- 
tor £ p can be expressed by £ p = 10s p r/~ 1 (0.1/e e ), e e is the 
fraction of the shock internal energy that goes into non- 
thermal electrons, r\ e is the radiative efficiency of these 
electrons, and <^ p is the proton acceleration efficiency. 
Modeling of the afterglows of gamma-ray bursts gives 
a typical value e e ~ 0.1 for relativistic shocks[12], so the 
typical value of £ p would be ~ 10. Thus the differential 



, dh 



£,pREx 



-3,-1 



ffi(^p,maa;/ £p,min) 



6xl0 4i £ Pi iergMpc"V 



(3) 

where e p ^ max and £p^ m in are, respectively, the maximum 
and minimum energy of accelerated protons and we have 
used ln(e P) max /e p , m m) = 10 in the last step. For £ p of the 
order ~ 10, this rate is much smaller than the required 
energy generation rate of cosmic rays per energy decade 
from 0.5 — 20 x 10 44 ergMpc _3 yr^ 1 deduced by different 
authors [l^. Note, however, that it was only the presence 
of short, powerful bursts early on that alerted us to its 
presence, so we can not exclude the possibility that many 
other similar, but rather less variable, events could be 
undetected. 

Pion production. - Now we consider the neutrino 
emission produced by these protons interacting with the 
soft photons in the sources. Protons lose ~ 20% of their 
energy at each pj interaction, dominated by the A res- 
onance. Approximately half of the pions are charged 
and decay into high energy neutrinos 7r + — > fi + + — > 
e + + v e + + v^, with the energy distributed roughly 
equally among the decay products. The fraction of en- 
ergy lost by protons to pions is f pl = R/(Tct ' ). De- 
noting by n(e 7 )c?e 7 the number density of photons in the 
energy range e 7 to e 7 + cie 7 , the cooling time of protons 
in the shock comoving frame for pj process is given by 



t> = 



—2 / dea{e)K(e)e / dxx' 

ZL pJeth Je/2T P 



2 n(x), (4) 



where T p = e p /Tm p c 2 , a and K are respectively the 
cross section and the inelasticity for p^f process [Til ]. 
The spectrum from infrared to X-ray frequencies of Sw 
J1644+57 can be approximately described by a broken 
power-law with n(e 7 ) = n{,(e 1 /ei y )~ a for e 7 < e& and 
n(e 7 ) = n;,(e 7 /eh) _ P for e 7 > et, where ~ IKeV, 
a ~ 2/3 and /3 ~ 2[4|. Approximating the integral by 
the contribution from the resonance we obtain 



f pi (£p) ~ 0.35- 



^X,47.5 



(ep/ept) 13 1 (e p < e pb ) 

Fffu^eb.KoV i ( £ p/ £ pb) a (s p > S p b) 



(5) 

where e pb = 0.15GeV 2 r 2 /e b = 1.5 x 10 16 r^e{; RcV eV is 
the proton break energy. To include the effect of multi- 
pion production and high inelasticity (which increases 
from ~ 0.2 at energies not far above the threshold to 
~ 0.5 — 0.6 at energies where multi-pion production dom- 
inates) at high energies [l5l[l6j. the above estimate of f pl 
should be multiplied by factor of ~ 2. As the neutrino 
energy is ~ 5% of the proton energy, the neutrino flux 
will peak at e vh ~ 7.5 x 10 14 r 2 ej; iReV eV. 

In the modeling of the SED of Sw J1644+57, upper 
limits from Fermi and VERITAS require T < 20 in the 
X-ray emitting region [4]. It is useful to express / P7 as 
a function of the pair production optical depth t 77 . The 
optical depth for pair production of a photon of energy Sh 

r 77 (eh) = = yij r(m e cV ' where l ii is the mean 
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free path. The fraction of energy lost by protons to pions 
is fpry ~ Yji ra pi^pii £ p/ e pb) a ^ 1 for protons with energy 
£ P > £ P b- Thus, we have[i~7| 

M e p ) * 0.5. 77 (100MeV) ( jj^^ ) Q£f ' • 

(6) 

Modeling of the SED of Sw J1644+57 requires 
r 77 (100MeV) > 1 [4], so we conclude that a significant 
fraction (> 50%) of the energy of protons accelerated 
to energies larger than the break energy, e p b ~ 10 16 eV, 
would be lost to pion production. 

Meson cooling. — The neutrino production efficiency 
will be reduced if the secondary mesons suffer from 
cooling before decaying to neutrinos and other prod- 
ucts [18j |. The pions and muons suffer from radiative 
cooling due to both synchrotron emission and inverse- 
Compton emission. The total radiative cooling time 
for pions is t'„ rad = Zm%£ l[te T m\<!jJ' B (\ + /i C )] ~ 
2 x 10 6 (e' 7r /lTeV)- 1 e B 1 _ 1 i u ; 1 48 i?? 4 . 5 r^ s, where U' B is the 
energy density of the magnetic filed in the shock re- 
gion and fic < 1 is the correction factor accounting 
for the inverse-Compton loss. By comparing this cool- 
ing time t'^ rad with the lifetime of pions = 7^r = 
1.9 x 10~ 4 (e^./lTeV) s in the shock comoving frame, 
where 7 W and r are the pion Lorentz factor and proper 
lifetime, one can find a critical energy (in the observer 
frame) for pions, above which the effect of radiative cool- 
ing starts to suppress the neutrino flux, 

e,,rad = 6.3 x 10 5 eB^iL~$rft t))2 TeV. (7) 

Similarly, by comparing the radiative cooling time t' rad , 
with the lifetime of muons r' , one can obtain a criti- 
cal energy for muons, above which the effect of radiative 
cooling starts to suppress the anti-muon neutrino flux, 

e^rad = 3-2 x 10 4 e-^ 2 1 £;, 1 4 / 9 rf^, 2 TeV. (8) 

The above estimates lead us to conclude that the neu- 
trino flux below ~ 10 16 eV is not affected by the meson 
cooling for typical parameters of TDFs. At higher ener- 
gies, however, pion cooling and muon cooling will sup- 
press the neutrino flux by a factor approximately given, 
respectively, by 

Gr = Mm{4„,/r;, 1}, C M = Min^^/r;, 1}. (9) 

The spectrum and flux of the neutrino flare. — The 
fluence spectrum of the muon neutrino {vp + Vp) emission 
from one TDF is 

£ 2<j, _ ,.2 dn P /p7Ctt(1 + (p) _ E X Cp/p 7 C7r(l + £ M ) 

(10) 

where e 2 p ^ = E p /[4wd 2 L ki(ep >max /e Pim i n )] is the differ- 
ential proton fluence produced by one TDF of total en- 
ergy E p = £ P -Ex in protons and the factor 1/8 repre- 
sents that the neutrinos produced by pion decay carry 
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FIG. 1: The expected muon neutrino (vp+Vp) spectra from 
Sw J1644+57 for different jet parameters. In all lines, Sx = 
7 x 10~ 4 ergcm- 2 , L x = 3 x 10 47 ergs _1 , t„ = 100s, e b = IKeV 
and €b =0.1 are used. 

one-eighth of the energy lost by protons to pion pro- 
duction, since charged and neutral pions are produced 
with roughly equal probability and muon neutrinos carry 
roughly one- fourth of the pion energy in pion decay. Fig.l 
shows the expected muon neutrino fluence spectra from 
Sw J1644+57, obtained by using the Lorentzian form for 
the photopion production cross section at the resonance 
peak plus a component contributed by multi-pion pro- 
duction at higher energies [2(| in calculating t' pl with 
Eq.(4). The initial rise in the spectrum at low energies is 
caused by the increasing pion production efficiency with 
energy, while the mild steepening and sharp steepening 
seen at higher energies are caused by muon cooling and 
pion cooling respectively. 

Now we estimate the number of neutrinos that 
can be detected from one TDF event similar to Sw 
J1644+57. The detection efficiency in water or ice 
of ultra-relativistic upward-going muon neutrinos with 
energies e v is P„ M ~ 7 x 10~ 5 (e„/10 4 - 5 GeV) K , where 
k = 1.35 for e v < 10 45 GeV, and n = 0.55 for e v > 
10 4 5 GeV[2l[. For neutrino fluence spectrum parameter- 
ized by £^$1,, the number of Vp and Vp above a certain 
energy e v $ detected by a km 3 -scale neutrino detector, 
such as IceCube, with area A = 10 10 ^4iocm 2 is 

N v (> £„ = 10 4 5 GeV) = / e ~ * v P»nAde v 

~ 2 Cp,1^10/p 7 ( £ P&) ( 3 xl0^crg) (lJfc) ' 

where we have used \n(e Pimax /e p<min ) = 10, ( w ~ 1, (p ~ 
1, a — 2/3 and = 2 in the last step. As / P7 > 0.5, 
we expect > 1 neutrinos above 30 TeV could be detected 
from Sw 1644+57 by Km 3 -scale detectors for £ p = 10. 
A careful calculation of the number of neutrinos above 
a certain energy as a function of the neutrino energy is 
shown in Fig. 2. A lower bulk Lorentz factor or a higher 
wind luminosity (i.e. a larger £ p ) leads to a larger number 
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FIG. 2: The expected number of neutrinos (f M +z> M ) above a 
certain energy detected from Sw J1644+57 by Km 3 -scale neu- 
trino detectors such as ICecube. The thin solid line represent 
the background atmospheric neutrinos. The same parameters 
as in Fig.l are used. 



of neutrinos that can be detected. If a similar event to 
Sw J1644+57 occurs at a closer distance (e.g. at z = 0.2) 
in future, more neutrinos would be detected as well. 

Neutrino detection from TDFs can be assured only if 
the number of background counts is smaller than one. 
The number of atmospheric neutrinos above 30 TeV ex- 
pected in the direction of the source during the flare pe- 
riod is 



N atm {> 10 
~3x 10~ 3 



4.5 



GeV) = J 



'AT' 
. 10 B s. 



10 45 GeV 
2 



de v Jdnjdt 



pair, 



(12) 

where AT is characteristic duration of the TDF, 8 ~ 
0.5° — 0.6° is the angular resolution of the neutrino 
detector at 30TeV-PeV[22|, Ff m is the cosmic-ray in- 
duced atmospheric neutrino background flux. We fit 
the atmospheric neutrino data measured by Icecube [23[ 
with a single power-law function, which gives F£ tm ~ 
4.7 x lQ- 8 ergcm- 2 s" 1 sr- 1 (e I ,/lTeV)" /3 with /3 ~ 1.44 in 
the energy range of 0.1-400 TeV. Since the number of at- 
mospheric neutrinos above 30 TeV expected in the direc- 



tion of the source during the flare period is much smaller 
than one, a detection of two neutrinos at energies above 
30 TeV from TDF sources will be highly significant. 

Discussions. - Neutrino emission has been pre- 
dicted to be produced by relativistic jets in gamma- 
ray bursts fill. 171 [24|, active galactic nuclei [HI and 
microquasars 



26j. Observations of Sw 1644+57 suggest 



that powerful jets are formed in TDFs, which have larger 
fluence in X/7-rays than the brightest GRBs and have 
higher X-ray luminosity than AGNs. There are three 
factors that are favorable for bright neutrino emission 
produced in such TDFs: 1) large fluence in the X-ray 
emission, which suggests large fluence in accelerated pro- 
tons; 2) high pion production efficiency as implied by the 
high opacity of high-energy gamma-rays, which leads to a 
high fraction of the proton energy lost into secondary pi- 
ons; 3) insignificant radiative cooling of secondary pions 
and muons, which leads to an almost flat neutrino spec- 
trum up to ~ 10 16 eV. The main uncertainty lies in the 
ratio between the energy density of protons and the en- 
ergy density in X-rays. In the magnetic field-dominated 
jet model, the proton energy density is subdominant, so 
the neutrino flux would be low, whereas in the matter- 
dominated jet model, we expect one to several neutrinos 
detectable by Km 3 -scale neutrino detectors from TDFs 
similar to Sw 1644+57. The Swift BAT, with a field of 
view of 4-7r/7 sr, has detected one TDF in 7 years, so the 
all-sky rate of TDFs would be one event in every one year. 
For neutrino detectors such as Icecube that has a 2-7T sr 
field of view, we expect one TDF event in the field of 
view of Icecube every two years, if the electro-magnetic 
counterparts can be identified. Therefore neutrino ob- 
servations provide a promising approach to diagnose the 
composition of the jets resulted from tidal disruption of 
stars by massive black holes in the galaxy center. 
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